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Preface to the Second Edition 

, In addition to the considerable interest of the results proved in this article, and the even greater 

interest of the methodology, this article is a landmark case in scholarly publishing. After the first 
version of this article was outright rejected by an anonymous referee of the Proceedings of the 
American Mathematical Society, because too many details were left to the reader (he or she didn't 
give us a chance to write a new version with more details), we decided to solicit nine non-anonymous 
CN ■ reports from world-class experts, assigning them specific parts. The division of labor (quite a few of 

them did over and above of what we asked them to, and refereed everything), and the full reports, 



> 

o 



. on the first edition, can be gotten from the webpage of this article: 

http://www.math.rutgers.edu/~zeilberg/mamarim/mamarimhtml/invmaj .html , 

o 

Q \ already mentioned in footnote 1. With a few exceptions of stylistic suggestions that we preferred 

not to adopt, we incorporated all of their (excellent!) suggestions. The first edition is also availabe 
there, for the record. 

X. 

H In order to be faithful to the original version, we have included in square brackets, and smaller 

font, the extra explanations demanded by the referees. We thank them profusely, and we now 
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believe that the formal correctness and clarity far exceeds %99.99 of the articles published in 
(anonymous!) "peer" -reviewed mathematical journals. Such journals, even electronic ones, will 
soon become obsolete, together with their pompous "editors" and anonymous referees. Instead, 
the present model of author (s)-appointed refereeing and self-publishing, in the authors' personal 
websites and the arxiv, with all the referee reports made public, and the referees acknowledged and 
given explicit recognition for their trouble, would become the norm, possibly with some tweaking. 
The present article is exclusively published in the Personal Journal of Shalosh B. Ekhad and Doron 
Zeilberger and arxiv . org. 

Human Statistics 

Human statistics are numerical attributes defined on humans, for example, longevity, height, weight, 
IQ, and it is well-known, at least empirically, that these are, each separately, asymptotically normal, 
which means that if you draw a histogram with the statistical data, it would look like a bell-curve. It 
is also true that they are usually joint-asymptotically-normal, but usually not independently so. But 
if you compute empirically the correlation matrix, you would get, asymptotically (i.e. for "large" 
populations) that they are close to being distributed according to a multivariate (generalized) 
Gaussian exp{—Q{xi, X2, ■■■)) with Q{xi,X2, ■ ■ ■) a certain quadratic form that can be deduced 
from the correlation matrix. 

Permutation Statistics 

Let our population be the set of permutations of {1, 2, . . . , n}. They too, can be assigned numer- 
ical attributes, and the great classical combinatorialist Dominique Foata (who got his Doctorat de 
troisieme cycle in statistics!) coined the term permutation statistics for them. 

The most important permutation statistic is the number of inversions, inv{'K), that counts the 
number of pairs 1 < i < j < n such that TTi > ttj (and ranges from to n(n — l)/2). For example, 
inf (314625) = 5, corresponding to the set of pairs {[1, 2], [1, 5], [3, 5], [4, 5], [4, 6]}. It features in 
the definition of the determinant, and Netto proved that the probability generating function (the 
polynomial in q such that its coefficient of g* is the probability that a uniformly-at-random n- 
permutation has i inversions) is given by 

(i)(i+g)(i+g+g^)---(i+g+g^+...+g"-^) _ nr=i(i-g') 

n! n!(l-g)" 

The second most important permutation statistic is the major index, maj{'n), that is the sum of 
the places i, where tt^ > TTj+i. For example, maj (314625) = 1 -|- 4 = 5, because at z = 1 and i = A 
we have descents. Major Percy Alexander MacMahon [M] famously proved that the probability 
generating function for the major index is also given by that very same formula. In other words 
the permutation statistics inv and maj are equidistributed. Dominique Foata [Fo] gave a lovely 
seminal bijective proof that proved the stronger statement that inv and maj are equi-distributed 
also when restricted to permutations ending at a given integer. 

William Feller ([Fe], 3rd ed., p. 257) proved that the number of inversions (and hence also the major 



2 



index) is asymptotically normal in the following sense. Feller easily computed the expectation, 



E[inv] = rUn = n{n — l)/4 

and the variance, 



2 _ 2n^ + - 5n 



If we denote by X„ the centralized and normalized random variable 

_ inv - rUn 
— , 



then Xn — 5- A^, as n — oo, in distribution, where N is the Gaussian distribution whose probability 
density function is e~^ /^/\/27r. 

A computer- generated proof, that gives much more detail about the rate of convergence to M, can 
be obtained using Zeilberger's Maple package 

http : //www .math. rutgers . edu/~zeilberg/tokhiiiot/AsymptoticMoments that accompanies the 
article [Z]. 

So both inv and maj are individually asymptotically normal, but what about their interaction^ 
In this article, we prove that they are asymptotically joint-independently-normal. In other words, 
defining, 

inv{n) -rUn , , . . maj{7r) - rrin 
A„(7rj := , r„(7rj := - 



we have that 

Pr{Xr, <s , Yn<t)^ ^ [ [ e-^'/2-yV2 ^y^^ n ^ oo. 

J —oo J —oo 



A Brief History 

It all started when the great Swedish probabilist Svante Janson (member of the Swedish Academy of 
Science, that awards the Nobel prizes) asked Donald Knuth (one of the greatest computer scientists 
of all time, winner of the Turing and Kyoto prizes, among many other honors) about the asymptotic 
covariance of inv and maj. Neither of these luminaries knew the answer, so Don Knuth asked one 
of us (DZ). DZ didn't know the answer either, so he asked his beloved servant, Shalosh B. Ekhad, 
who immediately ([E]) produced, not just the asymptotics, but the exact answerl It turned out 
to be n{n — l)/8. In particular, the correlation coefficient, Cov{inv,maj)/(T^ = 3„3^3^2_5„ = 

72 

^ + 0(l/n^) tends to zero as n goes to infinity. It followed that in the long-run, inv and maj are 
practically uncorrelated. 

But there are lots of pairs of random variables that are uncorrelated yet not independent. A 
convenient way to prove that and are asymptotically independent (we already know that 
they are both normal) is to use the method of moments, and to prove that the mixed moments 

Mr,s{n) := E[X:,Y^] , 
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tend to the mixed moments of A/" x M, as n — >■ oo . In other words, for r, s > 1: 



hm M2r,2s{n) 



(2r)! (2g)! 



hm M2r-i,2s{n) = , 
hm M2r,2s-i{n) = , 

n—^oo 

hm M2r_i,2s-i(n) = 



(EE) 

(OE) 
(EO) 
{OO) 



Ekhad's brihiant approach to the Janson-Knuth question merely settled the case r = l,s = 1 of 
(OO). Of course, because of symmetry {OE) and (EO) are trivially true (before taking the limits!, 
i.e. M2r,2s-iin) = and M2r-i,2s{n) = 0). 

[indeed if com ([vTi, 7r„]) := [n + 1 — TTi, . . . , n + 1 — TTjj] then trivially maj(7r) +maj(com(7r)) =n{n — 
l)/2 and inv{TT) + inv{com{TT)) = n{n - l)/2. So Pr({X„ = = j}) = Pr({X„ = -i,Yn = -j}) 
for all i, j, and it follows that 

M2r,2s-i{n) = ^iV*-^Pr({X„ = =i}) = ^iV*-iPr({X„ = -i,y„ = -j}) = 

Y.i-if'-i-jf'-^PriiX^ = i,Y^=j}) = -Y,{if^{jf^-^Pr{{X^ = i,Y^ = j}) = -M2r,2s-i{n) 

SO -M2r,2s— 1(^) equals to its negative, and so must vanish. The proof that M2r—l,2s{l^) = is similar.] 

One natural approach would be to extend Ekhad's brilliant derivation of Mi_i(n) to the general 
case, and try to derive closed-form expressions for M^^s{'n) for larger r and s. Since Ekhad's proof 
[E] is so brief, we can cite it here in full. 

"Svante Janson asked Don Knuth, who asked mc, about the covariance of inv and maj. The 
answer is (2) /4. To prove it, I asked Shalosh to compute the average of the quantity {inv{7r) — 
E{inv)){maj{TT) — E{maj)) over all permutations of a given length n, and it gave me, for n = 
1,2,3,4,5, the values 0,1/4,3/4,3/2,5/2, respectively. Since we know a priori^ that this is a 
polynomial of degree < 4, this must be it! □" . 

Obviously this brute-brute-force approach would be hopeless for deriving polynomial expressions 
for the moments Mr, sin) for larger r and s. As we will soon see, the degree of the polynomial 
M2r,2s{n) is 3(r + s), so for example, in order to (rigorously) guess Mio,io(w), we would need 31 

This is the old trick to compute moments of combinatorial 'statistics', described nicely in [GKP], section 8.2, by 
changing the order of summation. It applies equally well to covariance. Rather than actually carrying out the gory 
details, we observe that this is always a polynomial whose degree is trivial to bound. [Added in 2nd edition: the 
referees didn't find this obvious and asked for an explanation. See the bottom of page 8 and the top of page 9 in the 
present article.]. 
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data points, and we would have to ask our computers to examine more than 31! > 0.822 • 10^^ 
permutations. 



However, an inspired, still "empirical" (yet fully rigorous) "brute-force" approach does work. The 
first step would be to have a more efficient way to compute the moments Mr,s(«), for specific n 
and specific r and s. We will do it by first designing an efficient way to generate the probability 
generating function, let's call it G{n){p, q), for the pair of statistics (inv, maj). There are beautiful 
closed-form expressions for G{n){p,l) and G{n){l,q) (the same one, actually, due to Netto and 
MacMahon, given in page 2), but no such closed- form expression seems to exist for the bi-variate 
generating function, so the best that we can hope for is to find a recurrence scheme. 

A Combinatorial Interlude 

Let us forget about probability for a few moments, and focus on a fast algorithm for computing 

li-(n)(p,g) := ^ , 

for n up to, say, n = 50. 

[ Referee Dan Romik believe that we should mention, at this point, the "explicit" formula of Roselle[R] (mentioned 
by Knuth[K]) in terms of a certain infinite double product for the (/-exponential generating function of H(n){p, §). 
Romik believes that this may lead to an alternative proof, that would even imply a stronger result (a local limit 
law). We strongly doubt this, and DZ is hereby offering $300 for the first person to supply such a proof, whose 
length should not exceed the length of this article. 

Referee Christian Krattenthaler believes that we should also mention the beautiful extension of Roselle's result, by 
Adriano Garsia and Ira Gessel [GG], handling more permutation statistics. ] 

Define the weight of a permutation vr to be p^^vMq"^o-iM _ Suppose that tt G Sn ends with i, so we 
can write n = n'i, where vr' is a permutation of {1, . . . , i — 1, i -|- 1, . . . n}. 

When you chop off i from tt to form n' you always lose n — i inversions (that is, tt' has n — i fewer 
inversions than vr). The major index, however, decreases by n — 1 if the last letter of vr', let's call 
it j, is larger than i. 11 j < i the major index does not change at all. So writing vr = 7r"ji, we have 

inv{n"ji) = inv [it" j) +n — i , 

maj{TT"j), \{j<i; 
maj{'K" j) + n — 1, \i j > i. 



maj{iT"ji) 
Combining, we have 



weightijT" ji) 



p"- ^weight{n"j), if j < i ; 

p'^~^q'^~^weight{Tr" j), if j > i. 

So in order to compute H(n){p,q), we need to introduce the more general weight-enumerators of 
those permutations in Sn that end with an i. Let's call these F{n,i){p,q). In symbols: 
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It follows that (let's omit the arguments {p,q) from now on): 



i—1 n—1 
j=l j=i 

Note that, when we chop off the last entry, i, from tt = 7T"ji, 'k" j is a permutation of {!,...,? — 
1, i + 1, . . . , n}. We then "reduce" -k" j to a permutation of {1, . . . , n — 1} by diminishing all entries 
larger than i by 1. Hence the summation ranges from j = i\,oj = n—l rather than from j = i + 1 
to j = n. 

Replacing i by z + 1 in the above equation, we have: 

i n—1 

F{n,i + l)=p^-'-'Y.F{n-l,j)+p^-'-'q^-' J2 F{n-lJ) ■ 

j=l j=i+l 

Subtracting the former equation from p times the latter we get 

F{n, i) — pF{n, i + 1) = 

i—1 n—1 

i n — 1 

j=l j=i+l 

(i — 1 ^ \ I ''^"1 n — 1 

YF{n-l,j)-YFin-l,j) + p^-^' I ^F(n-l,i)- ^ F{n-l,j) 
3 = 1 j = l j \ j=i j=i + l 

= -p^'-'Fin - 1, i) + p^-'q''-^F{n - 1, i) = - l)F{n - 1, z) . 

Rearranging, we get: 

F{n,i) =pF{n,i + l) +p^-\q''-^ -l)F{n-l,i) for l<i<n . (RecF) 

We still need to specify F{n,n), and for this we do need the J2 symbol, namely we use Eq. {Fni) 
with i = n: 

n-1 

F{n, n) = Y F{n - 1, j) ■ (Fnn) 

The recurrence (RecF) together with the final condition (Fnn), and the trivial initial condition 
F{1, i) = 5i^i, enables us to efficiently compute F{n, i) for numeric (n, i), for {{n, < i < n < N} 
for any finite N (not too large, but not too small either: e.g., N = 100 is still feasible). In particular, 
we can compile a table of H{n){p,q) = F{n + l,n + l)(p, g), (the generating function for all n- 
permutations) for n < N — 1. 
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A crash course in mult i variable enumerative probability 

Suppose that you have a finite set of objects S and several statistics fi{s),. . . ,fr{s). The multi- 
variable generating function (weight-enumerator under the weig ) is defined to be: 

\" J^(') . . . rf-('^ 
ses 

Suppose that you pick an element s G S* uniformly at random and you want the multivariable 
generating function such that the coefficient of x"^ • • • x"'' would give you the probability that 
= oi, • • • , fris) = ar- It is giveu by: 

' ' ses 

The expectations, /i , . . . , are simply 

/. = (^p)(i.....i) . 

The centralized probability generating function is 

p/ \ -P(xi , . . . , Xr) 

Jr^i^Xi, ■ ■ ■ , Xr ) = ^ 

T-'-^ . . . ^■'^ 

The mixed moments (about the mean) of the statistics /i(s), . . . , fr{s) are defined by 

Mom[a,, ...,ar] = r^ Y.^fi{s) - hT' ' ' ' ifAs) - Uf^ . 
I-^I ses 

Often it is more convenient to consider the mixed factorial moments, using the combinatorial 
"powers" z^'^^ := z{z — 1) • • • (z — r + 1), better known as the falling-factorials. The mixed factorial 
moments are defined analogously by 

FM{a,, . . . , a,) = ^ " /i)^"'^ " " " (/-(^) " Z-)^""^ " 

Once you know the FM's for all ai,...,ar < M, you can easily figure out the Mom's (for 
ai,. . . ,ar < M), using Stirling numbers of the second kind (see [Z]). It is well-known and easy 
to see that one can just as well use the method of factorial moments in order to prove asymptotic 
independence. In other words, it would suffice to prove the analogs of (EE), (OE), (EO), (OO) 
with the moments E[X^Y'^] replaced by the factorial moments E[X^'^''Y^^^. 

The best way to compute FM{a\, ... , Qr) for all < ai, . . . , a,. < M, for some fixed positive integer 
M, is via the Taylor expansion of P{xi, . . . , Xr) around (xi, . . . ,Xr) = (1, . . . , 1), or equivalently, 
the Maclaurin expansion ofP(l + xi,...,l-|-Xr) 

ai,...,ay.>0 
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Back to inv-maj 



In the present approach, we need to put up with the more general discrete function F{n,i){p,q) 
even though ultimately we are only interested in H{n)(p, q) = F{n + l,n + l){p, q)- We will prove 
the stronger statement that even if you restrict attention to those (n — 1)! permutations that end 
with a specific i, the pair (inv, maj) is still asymptotically-joint-independently normal. 

Since the averages of both inv and maj over the permutations that end in i is n — z + (n — l)(n — 2)/4 
[for inv it is obvious, the last entry i contributes n — i to the number of inversions, and removing the last entry yields 
an n — 1-permutation, and for maj this follows from Foata's bijection mentioned above that maps inv to maj 
preserving the last entry], the centralized probability generating function corresponding to F{n, i){p, q) 
is: 

.w X _ F{n,i){p, q) 

^{n, t){p, q) .- _ j^,^^^^„_i+(„_i)(„_2)/4 ■ 

The recurrence {RecF) becomes 

G{n, i) = lG(n, i + l) + ''^^S^l 2 1) " ^' ' (^^^^) 
and the final condition becomes 

^ n — l 

G(n,n) = — -5^(Mr/2-^G(n-l,i) . (Gnn) 

We also need the obvious initial condition G{l,i) = (5j^i. 

Guessing Polynomial Expressions for the Factorial Moments 

Equipped with these efficient recurrences, our computer computes G{n,i){p,q) for many values of 
n and i. Then for each of these it computes the (r, s) (mixed) factorial moments FM{r, s)(n, i), for 
many small numeric (r, s) by computing the (initial terms of the) Maclaurin series for G(n,i)(l + 
p,l + q). We then fix numeric values of (r, s) and use polynomial interpolation to guess explicit 
polynomial expressions for FM[r, s){n,i) as polynomials in (n, i) for that pair (r, s). The process 
is repeated for all pairs (r, s) for which < r, s < M, for some pre-determined specific positive 
integer M. Wc note that it is obvious, both from the combinatorics and from the recurrences, that 
the FM(r, s){n, i) are always polynomials in (n, i), for any fixed numeric r and s. 

[As we have already mentioned in footnote 2, most of the referees didn't find this obvious. The proof via the 
recurrences {RecG') and (Gnn') to be derived in page 10 is by induction on (r, s) and the fact that the indefinite 
sum of a polynomial (in this case with respect to i) is yet another polynomial, and the "operator" on the right of 
(Gnn') is polynomial-preserving. 

Let's sketch the combinatorial proof that the mixed moments Af^. ^(n) are polynomials in n. The combinatorial 
proof that FM(r,s){n,i) are polynomials in both n and i is similar. Write inv{Tr) and maj{lT) as a sum of 
"atomic" events, e.g. for inv the sum of xij^j > ) over all pairs of integers (i, j) satisfying 1 i < j n. Here 
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X{S) = 1 if is true and x(S') = if it is false, maj can be similarly expressed as a sum of < j AND TTj > 
TTj-i-i). The sum of inv(Tr)^maj(Tr)^ over all permutations vr € Sn can be expressed as a multi-sum with the 
outer sum ranging over and the inner sums with 2(r + s) sigma signs. Now do discrete Fubini! Bring the 
formerly outer-sigma, over S^, all the Way inside past all the other 2(r + s) sigma signs. Each individual sigma 
sign involves two indices, and collectively these 2(r + s) sigma signs involve < 2(r + s) indices, corresponding 
to locations in a generic 7T,-permutation TT, some of whom may coincide. Let's call them ii < ... < ife ( where 
A; < 2r + 2s). These can be placed in lots of possible intertwining ways, and so can the values of TT in those 
places. There are still finitely many scenarios. (Formally, one gets a Cartesian product of two partially ordered sets, 
one for the domain and one for the range, each of which has finitely many linear extensions. This is reminiscent of 
Richard Stanley's theory of P-partitions). For each particular such scenario (linear extension of the domain-poset 
and the range-poset) , there are ways to choose the participant indices, and ways to choose their occupants 
(i.e. the values of TT there) and the remaining n — k entries can, of course, be arranged in {n — k)\ ways, yielding 
(fc) ~ ways. Dividing by n\ gives (^) {n — k)\/n\ = (^)/A;!, a polynomial in n of degree k. Now the 
whole thing is a sum of finitely many such {D /k\ for < A; < 2r + 2s, and since a finite sum of polynomials is 
still a polynomial, we are done! Now isn't that ohviousV. ] 

It would have been nice if we could guess closed-form expressions for FM{r, s){n,i) for symbolic 
(r, s), but no such closed-form exists as far as we know, and besides it is too much to ask for 
and more than we need. To prove asymptotic normality we only need the leading terms. Viewing 
the leading terms, our beloved computer easily conjectures the following expressions. For integers 
r > 0, s > we have: 

FM{2r,2s){n,i) = Vn"Vi ^ n^""^^" + {lower - total - degree - terms - in - {n,i)) . 



2^r[ 2«s! V36, 

[Note that the coefficients of n^'^~^^^~^i, n^^~^^^~^i^ etc. are all zero, hence they don't show up!]. 
For integers r > 0, s > 1 we have: 

FM{2r,2s - l)(n,i) = n^''+^'-^[-{s - l)n^ - 6rnH + ISrm^ - 12ri^] 

+ {lower — total — degree — terms — in — {n, i)) 

For integers r > 1, s > we have: 

r+s — l 



(2r)! (2s)! f 1 
~¥r\lFs 

For integers r > 1, s > 1 we have: 

(2r)! {2i 

2^2^s! V36y 2 
Nice conjectures but what about proofs? 



FM(2r— 1, 2s)(n, z) = — -^^J-^-^^p^ ( ^ ) {i — l)n^'^^^^ ^ + {lowei — total— degree— terms— in— {n,i)) . 



(2rV (2s)' / 1 \ '^+*~^ 9 

FM(2r-l,2s-l)(n,i) = Vn"V^ ^ -n^'"'^^~^ {n-2if + {lower-total-degree-terms-in-{n,i)) 
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While we prefer the empirical approach of guessing, an alternative approach to finding many 
FM(r, s){n,iys, that is also necessary in order to rigorously prove our conjectures, is to first use 
(RecG) and (Gnn). Write G{n,i){l + p,l + q) as an infinite generic Taylor series around (0,0), 
and write down the implied infinite- order recurrences expressing FM{r,s) in terms of FM{r',s') 
with r' + s' < r + s. 

[The infinite-order recurrence for the FM(r, s)(n,i), obtained from {FiecG) is gotten by expanding 1/(1 + q) 
and 

{\+pY-i{l + qY-i _ 1 

((l+p)(l + g))"/2 

as Maclaurin series in {p,q), using the binomial theorem and manipulations on formal power series (that Maple 
does automatically to any desired order), and combining terms. Similarly, the implication of {Gnui) is obtained by 
expanding ((1 +p)(l + q))^^^~-' using the binomial theorem. Both of these tasks are accomplished by procedure 
MOP in the Maple package InvMaj . ] 

Note that in order to compute FM{r, s){n,i), for any specific, numeric r and s, we only need 
finitely many terms (actually rs — 1 of them) of the infinite-order recurrence, since eventually all 
the contributions will be zero. Of course, as we have already commented, there is no hope for 
finding a general expression for FM{2r,2s){n,i) , FM{2r,2s — l){n,i) , FM{2r — l,2s)(n,i) and 
FM{2r — 1,2s — l)(n, i), depending explicitly on r and s (i.e. symbolically in terms of r and s) 
as well as on n and i, but to prove, by induction on r, s, that the above leading terms are valid, 
all we need to do is to verify that the leading terms of the implied recurrences for the FM{r, s)'s 
(that we have just talked about) are consistent with the above explicit expressions. 

The implication of (RecG) is 

FM{r, s)(n, i) - FM(r, s)(n, i + = 

-sFM{r, s-l){n,i + l) + sFM{r, s-l){n-l,i) + — '-FM{r - 1, s - l)(n - 1, z) + 

{lower — order — terms) , {RecG') 

while the implication of {Gnn) is: 

n— 1 n— 1 

2(n-l)^' 
j=i ^ ' j=i 



FM{r,s){n,n) = ^ FM{r, s){n -l,j)- ^(2j - n)FM{r, s - l)(n - 1, j) 



n—1 n—1 



— — - ^(2j - n)FM{r - l,s)(n - l,j) + ——- Y,{2j - n)^FM{r - l,s - l){n - l,j) 

+{lower — ordei — terms) . {Gnn') 

[Referee Guoniu Han correctly commented that one must say something about the degree of the polynomials in (n, i) 
that reside in the "lower-order terms" that feature as coefficients in these infinite-order recurrences. It turns out the 
the coefficient of FM{r — r' , s — s'){n — 1, i) in {RecG') and the coefficient of FM{r — r',s — s'){n, i + 1) 



10 



(and similarly for [Gnn')) are polynomials in (n, i) of total degree < r' + s', thanks to the binomial theorem, and 
this would imply, by induction on (r, s) that not only are FM(2r,2s)(n,i) etc. polynomials, but their leading 
terms look as claimed above.] 

The next step is to spell out these two recurrences each into the four cases according to whether 
(r, s) is {even, even), {even, odd), {odd, even), and {odd, odd). 

Once you have these eight recurrences, for each and every one of them, you plug in the above 
conjectured expressions for the leading terms, and verify that up to the leading terms, things agree. 
At the end of the day, after dividing by 

(2r-2)! (2s -2)! 
(r- l)!2^-i (s- l)!2«-i ' 

this boils down to proving equalities among certain low-degree polynomials in (r, s) (namely the 
leading coefficients, in {n,i)), that in turn, reduces (since A = B iS A — B = 0) to proving that 
certain low-degree polynomials in (r, s) are identically zero. 

So in order to check that these low-degree polynomials in (r, s) are all identically zero, it is enough to 
check each of them for finitely (and not-too-many) numeric r, s. Typing Checkl (FM8m(n, i) ,n, i) ; 
and Check2(FM8m(n,i) ,n,i) ; in the Maple package InvMaj does exactly that, by checking that if 
you plug in the conjectured leading terms of FM{2r,2s){n,i), FM(2) — l,2s){n,i), FM{2r,2s — 
l){n,i) and FM{2r — 1,2s — l){n,i) and subtract the right sides from the left sides (for each of 
the eight cases) you get lower-order polynomials, in {n,i), for all 1 < r, ,s < 8. This proves all 
these claims (rigorously), with a vengeance! The (8/2)^ = 16 special cases are much more than is 
needed, since the relevant polynomials in (r, s) are easily seen to have degree < 2 so (2 + 1)^ = 9 
agreements would have sufficed. 

The Maple package InvMaj 

All the nitty-gritty calculations described above, that constitute a fully rigorous proof, may be 
found in the Maple package InvMaj accompanying this article. This package is available from the 
webpage of the present article: 

http : //www . math . rutgers . edu/~zeilberg/mamarim/mamariinhtml/invmaj . html , 

where the reader can also find some sample input and output. The direct url of the package is: 
http : //www .math . rutgers . edu/~zeilberg/tokhniot/InvMaj 

La Grande Finale 

The special cases r = 1, s = and r = 0, s = 1 give 

FM{2,0){n,i) = ^n^ + O{n^) , 

FM{0,2){n,i) = ^n^ + O{n^) . 
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So (recall that we are interested in the normalized mixed factorial moments) 

FM{2r,2s){n,i) (2r)! (2s)! 



FM{2,0){n,iyFM{0,2){n,iy 2^r\ 2«s 
FM{2r,2s - l)(n,i 



+ 0(l/n) , 



FM{2,0){n,iyFM{0,2){n,iy-'^/^ ^ ' 
FM{2r-l,2s){n,i) 



FM{2, 0)(n, i)^- V2i?M(0, 2)(n, i) 
FM{2r-l,2s-l){n,i) 



o(l/n) , 



™(2,0)(n,z)^-V2FM(0,2)(n,z)«-V2 W J • 

And we see that as n ^ oo these indeed converge to the mixed moments of the famous mixed 
moments of the bivariate independent normal distribution e~" ^^/{2ir) □. 

Encore: A more refined asymptotics for the (Normalized) Mixed Moments 

With more effort, we (or rather, our computer) can gucss-and-prove the following asymptotics for 
the case of interest (n + 1, n + 1), i.e. the asymptotic expressions for the centralized-and-normalized 
(genuine, not factorial) mixed-moments, let's call them a(r, s)(n), for the pair of random variables 
{inv,maj) acting on the set of permutations of length n. Indeed, according to S. B. Ekhad, we 
have: 

(2r)! (2s)! / 9(r2 + _ ^ _ i i 



a(2r,2s)(n) = ^^^^ 1 - ^ —-- ■ - + 0(^' 

V ' A y 2'-r! 2^s! V 25 n ^n^' 

N (2r)! (2s)! / 9 / 81,2 2. 243, , 1773\ 1 ^,1, 

Of course, by symmetry a{2r, 2s— l)(n) and a(2r— 1, 2s)(n) are identically (not just asymptotically!) 
zero. 
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